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I We investigate the single-point probability density function of the velocity in three- 

dimensional stationary and decaying homogeneous isotropic turbulence. To this end we 
apply the statistical framework of the Lundgren-Monin-Novikov hierarchy combined 
with conditional averaging, identifying the quantities that determine the shape of the 






i-rt probability density function. In this framework the conditional averages of the rate of 

energy dissipation, the velocity diffusion and the pressure gradient with respect to velocity 
play a key role. Direct numerical simulations of the Navier-Stokes equation are used to 
complement the theoretical results and assess deviations from Gaussianity. 



I 

^. 

O 

^ 1. Introduction 

>% 

(-H The spatio-temporal complexity of turbulent flows demands a statistical description, 

Q^ which can be formulated in terms of probability density functions (PDFs) of the fluctu- 

ating turbulent quantities such as velocity or vorticity. The Lundgrcn-Monin-Novikov 
(LMN) hierarchy ( |Lundgren|1967[|Monin|1967[|Novikov|1968[[Uhnich fc Lyubimov|1969 1 



J> provides a statistical framework, which, starting from the basic equations of motion (i.e. 

On the Navier-Stokes equation or the vorticity equation), yields the temporal evolution of 

^^ the multi-point PDFs of the velocity or vorticity. As typical for a statistical theory of 

^«. turbulence, this approach has to face the famous closure problem of turbulence. When 
deriving the evolution equation of the single-point PDF, this equation couples to the two- 

trr* point PDF. Now determining the evolution of the two-point PDF involves the three-point 

^^ PDF and so forth, ending up in an ever increasing number of evolution equations for the 

,_^ multi-point PDFs. This chain of equations may be truncated at a given level by introduc- 

IL" ing reasonable closure approximations, which express the {n + l)-point PDF in terms of 

. ^H the n-point PDF. However, these approximations have to be performed with great care as 

j^ a too crude approximation results in serious defects. An alternative approach is provided 

M by conditional averaging, which, for example, is taken as a starting point for modelling 

C^ 



in the realm of the well-known PDF methods by Pope (2000). Novikov (19931 used this 
method in the case of the vorticity equation to elucidate the statistical balance of vortex 
stretching and vorticity diffusion. When combined with the LMN hierarchy, a truncation 
at a given level is possible by introducing the conditional averages of the terms appearing 
in the equations of motion as unknown functions. These functions then may be modelled 
or estimated from experiments or direct numerical simulations as recently exemplified in 
(Wilczek fc Friedrich|2009 1 in the case of the single-point vorticity PDF. The advantage 



of this approach lies in the possibility to interpret the unknown functions in terms of their 
physical meaning. For instance, it turns out that the single-point velocity PDF depends 
on the conditionally averaged rate of energy dissipation. 
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Ever since the experimental measurements by Townsend presented in (Batchelor|1953l 



it is an ongoing discussion, whether the single-point velocity PDF displays a Gaussian 
shape. Arguments in favour of this hypothesis often involve the central limit theorem, 
taking the velocity field as a spatial ensemble of independently fluctuating random vari- 
ables with finite variances. This argument turns out to be rather hand-waving as the 
assumption of statistical independence is obviously violated. Moreover, the same kind of 
argument should be applicable to the vorticity PDF. This PDF, however, is known to 
display extremely non-Gaussian tails, which are often associated with the existence of 
vortex structures in the flow. When applied with care, however, the central limit theorem 
together with proper assumptions on the energy spectrum has been shown to produce 



sub-Gaussian tails by Jimenez ( 1998 ) 



While the assumption of Gaussianity at least seems to be a good approximation, ex- 
perimental and numerical measurements indicate the possibility of sub-Gaussian tails 
([Vincent fc"Mencguzzi|1991[|Noullez et a/.|1997{|Gotoh et a/.|2002J. Theoretical support 



in this direction was, in the case of forced turbulence, given by Falkovich fc Lebedev 



(1997). The approach, however, depends on the type of forcing applied to the flow. This 
point recently was raised in ( iHosokawa^^2008 ) , where it was hypothesized that decaying 
turbulence might be profoundly different from forced turbulence. In this work, a closure 
for the LMN hierarchy was suggested leading to a Gaussian velocity PDF for decaying 
turbulence. Gaussian PDFs are also found as a result of the cross-independence hypoth- 



esis (Tatsumi & Yoshimura 2004) for both the single-point PDF and joint PDFs in the 
inertial range. In both works the pressure contributions are neglected on the level of the 
single-point PDF. 

The present work aims at clarifying these points. We make use of the statistical frame- 
work of the LMN hierarchy combined with conditional averaging to identify the terms 
which determine the shape and evolution of the single-point velocity PDF. This approach 
gives concise insights how the different terms interplay to determine the functional form 
of the PDF. The theory applies to forced as well as decaying turbulence. The results pre- 
sented in (Hosokawa 2008) can be obtained via a simple approximation, which neglects 
the statistical dependence of the pressure gradient and the energy dissipation on the 
velocity. Testing the theory against results from direct numerical simulation, however, 
indicates slightly non-Gaussian PDFs for both forced and decaying turbulence, which is 
related to the statistical correlations between the velocity and the different dynamical 
effects such as pressure gradient and energy dissipation. 

The remainder of this article is structured as follows. After reviewing the LMN hier- 
archy and conditional averaging, we discuss the implications of statistical symmetries, 
which help to drastically simplify the description. We then derive explicit formulae for 
the homogeneous and stationary velocity PDF. The conditional averages appearing in 
these formulae have to obey certain constraints which will be derived and used to moti- 
vate a simple closure approximation. We then turn to the numerical results of stationary 
and decaying turbulence examining the different quantities arising in the theoretical de- 
scription. To highlight the genuine properties of the velocity statistics, some comparisons 
with the single-point statistics of the vorticity are drawn before we conclude. 
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2. Theoretical Framework 

2.1. The Lundgren-Monin-Novikov Hierarchy 

As an introduction into the statistical framework used in this work, we start with re- 
visiting the derivation of the LMN hierarchy. For further details, we refer the reader to 
(|Lundgren"1967| |Monin|[l967| |Noviko^[l968l [Pope||2000 1 . 

The Navier-Stokes equation for an incompressible fluid reads as 

d 



dt 



u{x, t) + u{x, t) ■ Vuix, t) = -Vpix, t) + iyAu{x, t) + F{x, t). 



(2.1) 



Here u denotes the velocity field, p is the pressure, v denotes the kinematic viscosity and 
F denotes a large-scale forcing applied to the fluid to produce a statistically stationary 
flow. In the case of F{x,t) = 0, turbulence decays due to dissipation and statistical 
stationarity cannot be maintained. 

Let u{x,t) denote a realization of the velocity field and v denote the corresponding 
sample space variable. The fine-grained PDF of the velocity is then defined by f{v; x, t) = 
6{u{x,t) — v), from which the PDF can be obtained by ensemble-averaging. 



fiv;x,t) = {f{v;x,t)) = {6{u{x,t) ^ v)). 
Taking the derivative of the fine-grained PDF with respect to time yields 



d_ 
dt 



f{v;x,t) ^ -V„ 



du 
~dt 



{x,t)f{v;x,t) 



(2.2) 



(2.3) 



Together with an analogous calculation for the advective term we obtain 



d_ 
Ft 



/ + V- uf 






{-Vp + vAu + F) f 



(2.4) 



Here, the Navier-Stokes equation (2.11 has been used to replace the terms of the left- 
hand side with the terms of the right-hand side. Incompressibility (V • it = 0) was used 
to rearrange the advective term. In order to turn this into an equation for the PDF, 



we have to perform an ensemble average of (2.4). The averaging of the left-hand side is 
straightforward yielding dtf + v ■ ^ f; however, the terms on the right-hand side cannot 
be expressed in terms of / and v only. Now there are several options how to proceed. 
One possibility is to express these unclosed terms on the right-hand side with the help 
of the two-point PDF. We exemplify this for the pressure term. 

The two-point PDF of the velocity can be written as f2{vi,V2;xi,X2,t) = {S{u{xi,t)~ 
Vi) 6{u{x2,t) — V2))- Now the pressure term can be expressed in terms of the velocity 
according to 



((-V.,P)/) 



1 

An 

1 

An 



da;Jv.,^^ 



[u{x2,t) ■Va;2u{x2,t)] 



dX2 Va;i 



\X2 -Xi\ 
TT[{V^,Vljiu{x2,t)u^ix2,t))] ' 
\X2 -Xi\ 



S{u{xi,t) — vi) 



S{u{xi,t) ~ Vi) [ / dV26{u{x2,t) - V2] 
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All 



tTTj V \X2-X1\J 

(2.5) 



where we have assumed an infinite domain without further boundary conditions. Here, 
the identity J dw2 (5(ii(a;2, i) — V2) = 1, incompressibihty and the sifting property of the 
delta function, u{x2,t) 5{u{x2Tt) — V2) = V2 S{u{x2,t) — V2), have been used. This manip- 
ulation explicitly shows how the single-point statistics couples to the two-point statistics. 
The remaining unclosed terms can be treated in a similar manner; for more details, we 
refer the reader to (Lundgren||1967). Now determining the evolution equation for /2 ba- 



sically involves the same steps, leading to a coupling to /s. This eventually leads to the 
aforementioned never-ending hierarchy of evolution equations for the multi-point velocity 
distributions. 

Instead of establishing this chain of evolution equations, it is possible to truncate the 
hierarchy on a given level by introducing conditional averages according to e.g. 

((- Vp)/) = {{^Vp{x, t) ) S{u{x, t)-v))^{- Vp\v)f, (2.6) 

i.e. now the coupling to a higher-order PDF is replaced by introducing unknown functions 
in the form of conditionally averaged quantities, which may be modelled or obtained by 
experimental or numerical means. The conditional averages may be regarded as a measure 
for the correlation (in the sense of statistical dependence) of the different quantities on 
the single-point level. We will focus on this approach in the following. The temporal 
evolution of the velocity PDF then takes the general form 

^/ + t; • V/ = -V, • [( - Vp + lyAu + F\v)f] . (2.7) 

To rewrite the advective term, here again incompressibihty and the sifting property 
have been used. In this kinetic description for the evolution of the single-point velocity 
PDF, several unclosed terms appear: the conditionally averaged pressure gradient, the 
Laplacian of the velocity as well as the external forcing. We take this equation as a 
starting point for the following considerations and will now utilize statistical symmetries 
to simplify the description. 

2.2. Statistical Symmetries 
2.2.1. Hom,ogeneity 

For homogeneous flows, several simplifications arise. First, the advective term on the 



left-hand side of (2.7 1 vanishes as homogeneity means that / cannot depend on x. This 
simplifies the kinetic equation resulting in 

^^f = -V^-[{-Vp+,yAu + F\v)f]. (2.8) 

Note that in the case of decaying turbulence, the evolution is determined by the condi- 
tionally averaged pressure gradient and the dissipative term only. Homogeneity allows to 
recast the latter. Calculating the Laplacian of / and averaging yieldiljj 

d^ „ . 9/aV: \f+^f(^p^v)f. (2.9) 

dvi dvj \ dxk axk / 



dxl, dvi \ dx^ 



f Throughout the paper the Einstein summation convention is assumed. 
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This is a multi-dimensional version of the homogeneity relation introduced by |Ching| 
(1996). With this expression, the kinetic equation takes the form 



d_ 
dt 



f 



d_ 

dvi 



d 



da 



F, 



f 



_d__d_ 

dvi dvj 



dui duj 
dxk dxk 



f, 



(2.10) 



This 



which has also been derived for turbulence with mean flow by Pope (2000). 
troduces the tensor D^ = (i'{dui/dxk){duj/dxk)\v), sometimes termed conditionally 
averaged dissipation tensor, which will be of central interest in the following. 

Before we proceed to the simplifications that arise due to isotropy, we pause for a side 
remark as we would like to highlight the connection to moment equations which can be 



obtained directly from the PDF equation (2.101. The law of energy decay may serve as 



an example. In the decaying case, the only two terms in (2.10) governing the evolution 



of the velocity PDF are the conditionally averaged pressure gradient and the conditional 
dissipation tensor. While the former takes the form of a drift induced by the non-local 
pressure contributions, the latter term may be interpreted as a diffusive term, however 
with a negative sign. This sign can be explained in a physically sound way as for a 
diffusion process an initially localized concentration spreads over time, eventually being 
dispersed over a large domain. Regarding the decay of the velocity field of a turbulent 
fluid, the opposite takes place. An initially broad distribution of velocity contracts as the 
velocity field dies away. When the fluid has come to rest, the PDF is localized sharply in 
probability space, lim f{v;t) — S{v), expressing that we have probability one finding a 

vanishing velocity. 



d 



To proceed, we multiply (2.10) by v'^/2 and integrate over v, 

2 d 



at 



E, 



kin ' 



= / dv 



V 



2 dt 



f = -dv 



VkVk 



d 



-P 



v)f + 



d d 



dui duj 



dvi dvi \ dxh dxh 



v)f 



d_ 

dvi \ dxi 

■(2.11) 

The right-hand side may be integrated by parts. Assuming a sufficiently rapid decay of 
the PDF, one finds 



d 

dt' 



= / dv 



Vp\v)f 



dui dui 
dxk dxk 



= -{u. 



Vp, 



-> 



(2.12) 



The validity of the second equality will become clear in the next section (see (2.23)) 



The last equality comes from the fact that the pressure-related average vanishes due to 
homogeneity and incompressibility and that the rate of energy dissipation and squared 
vorticity (multiplied by v) have the same spatial and hence also ensemble average. Equa- 
tion (2.12) is what one expects and also what is long known as the evolution equation 



of the kinetic energy. It is not surprising that the correct equation for the energy fol- 
lows from the equations for the velocity PDF as those were derived directly from the 
Navier-Stokes equation without using any approximations. 

2.2.2. Isotropy 



The conditional averages appearing in (2.8) and (2.10) are vector- and symmetric 



tensor- valued functions of the velocity vector v. Statistical isotropy, which is here assumed 
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as invariance under rotatiorjf] imposes further constraints on these functions e.g. 

D{Kv) = RD(t>)R'^ R e S0(3), (2.13) 

where R is a rotation matrix. In three dimensions it can be shown that, because of 
these constraints, isotropic vector- and symmetric-tensor- valued functions hke (Vpv) 
and T){v) take the form ( |Robertson||1940| |Batchelor||1953[ ) 

a,{v)^a{v)'^, (2.14) 

V 

B,jiv) - ^l{v) 5,, + [\{v) - ii{v)] '^, (2.15) 

where a{v) = i) ■ a{v) is the projection of a{v) onto the unit vector in the direction of 
the velocity and A(w), n{v) are the eigenvalues of the matrix B. Here A is the eigenvalue 
of the eigenvector v^ whereas ^ is the eigenvalue of the (two-dimensional) eigenspace 
perpendicular to v. Note that because a, A and /i are isotropic scalar-valued functions, 
they depend only on the absolute value of the velocity v. The same goes for the probability 
density f{v), i.e. it is only a function of the absolute value v. This, however, may not be 
confused with the PDF of the absolute value of the velocity, in the following denoted by 
f{v). However, there is a simple relation between those two, 

f[v)^ATTv^f[v), (2.16) 

indicating that it suffices to determine the PDF of the absolute value in order to specify 
the PDF of the full vector. 

From the consideration about isotropic vector-valued functions, we find 

{-Vp\v) =Jl{v)v, H(i;) = (-M- Vp|w), (2.17) 

(i^AMJt;) = A(t;)t), A(v) = (i^it • Am|u) , (2.18) 

{F\v) ^^(v)v, $(w) = (m-F|w), (2.19) 

where ii and v denote the unit vectors in the direction of u and v. Note that it is sufficient 
to take the conditional average with respect to the magnitude of the velocity to obtain 
the functions H, A and $. This considerably simplifies the numerical estimation of the 
conditional averages. Let us now turn to the conditional dissipation tensor, which takes 
the form 

V>{v)^{vkk^\v), (2.20) 

where Aij = dui/dxj is the velocity gradient tensor. A may be decomposed into sym- 
metric and antisymmetric parts according to A = S -I- W, where S = ^(A -I- A""") and 
W = 5 (A — A"""). These two tensors characterize the local rate of stretching and the 
rate of rotation of the fluid. In the case of statistical isotropy, D is determined by its 
eigenvalues (see ( |2.15 )), which can be obtained using the relations 

Tr(D) = A(i>)-f 2/i(z;), (2.21) 

vD* = X{v). (2.22) 

The trace of D is determined by the conditional averages of the local rate of energy 
dissipation e = 2i/Tr(S^) and the squared vorticity w^: 

Tr(D) = i/(Tr(S2)-Tr(W2)|t;) ^ ^ {e + iyu^\v) . (2.23) 

f Invariance under reflections is not necessary for the results presented in this paper. 
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The second scalar quantity needed to determine D is 

vBv = {iyuAA^u\v) = {i^{A'^u)'^\v) . 



(2.24) 



This rather formally looking quantity has a simple physical interpretation. As A may 
be decomposed in symmetric and antisymmetric parts, we write A'^u = (S — W)u. The 
last term may also be written as Wii = ^u) x ii due to the relation Wij = —\eijk(jJk- 



Hence, the conditional average appearing in (2.24 1 involves the absolute value of the 



difference between the rate of stretching in the direction of the velocity vector and rate 
of rotation of the unit vector ii. Summing up, the conditional dissipation tensor D in 
isotropic turbulence has the form 



\iv) 



1 , 



[Xiv) 









^uf\v) 



(Kat«)» 



(2.25) 

(2.26) 
(2.27) 



These relations can now be used to simplify the structure of the PDF equations (12. 81 



(2.101. It turns out that by exploiting statistical isotropy, we have to deal with an effec- 



tively one-dimensional problem. Take for example 



_d_ 

dvi 



{vAui\v)}{v) 



dvi 47^1-^ 



Rv) 



1 d 



Hv)f{v), 



(2.28) 



Attv^ dv 

which makes clear that this term depends on v only. The same applies to terms involving 
D, a short calculation yields 



JLJL 
dvj dvi 



A,(^)/(t^) = 



1 



Attv^ 






X{v)f{v) 



d_2 

dv V 



tJ'iv)f{v) 



df' 



This leads to the isotropic form of (2. 

_a_ 

dv 
d_ 

dv 
d_ 

dv 



(2.10) 



= 



dV 



(n + A + $)/, 

v 



A- 



n + $ 



/ 



2^/' 



dv'^ 



A/. 



(2.29) 

(2.30) 
(2.31) 
(2.32) 



We note in passing that these equations can also be derived without the assumption of 
statistical isotropy and therefore correctly describe f{v) also in statistically anisotropic 
turbulence. However, in that case they are not equivalent to (2.8)-(2.10) and f{v) is 



not fully determined by f{v). Furthermore, the relations ( 2.17 1-( 2.19) and (2.25) are not 
valid in that case and do not allow a simple interpretation of the quantities H, A, $, A 
and /i. 

2.3. Homogeneous and Stationary PDFs 



Because now having reduced the problem to a one-dimensional, ( 2.31 ) is easily integrated 
yielding the velocity PDF in homogeneous turbulence. 



/>;i) 



N 
X{v,t) 



exp 



dv 



,A{ v\t) + ^ii{v\t) 
X{v',t) 



(2.33) 
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showing that it can be expressed as a function of the conditional averages involving the 
velocity diffusion, dissipation, enstrophy and the stretching and turning of the velocity 
vector. Here A/" denotes a normalization constant, which depends on uq. Note that by 
only imposing isotropy and homogeneity, the conditional averages will be a function of 
time. In the following, we will refer to this solution as the 'homogeneous solution'. 
In the case of forced turbulence, it is possible to maintain a statistically stationary 



flow. This immediately implies dtf — and from (2.30) follows 



n{v) + K{v) + ^{v) =0, 



(2.34) 



i.e. the pressure gradient term, the viscous term and the term stemming from the external 
forcing identically cancel. This also means that in the case of stationary, homogeneous 



and isotropic turbulence the probability current on the right-hand side of (2.8 1 identically 



vanishes. Note that this is not obvious and is due to the fact that the problem becomes 
effectively one-dimensional by virtue of isotropy. In the case of the turbulent vorticity, 



this balance was discussed in (Novikov||1993 Novikov & Dommermuth 



recently in (Wilczek & Friedrich 2009"). In the stationary case, (2.31) and (2.32) are 



1994) and more 



equivalent due to the conditional balance (2.34). The stationary solution of (2.32) reads 
as 

A/- r . , -nK) - ^{v') + |mK) 



/(«) 



exp 



dv' 



(2.35) 



and can be either obtained by solving (2.32) or from ( |2.33| ) by utilizing the conditional 
balance and omitting the time dependence. We emphasize here that introducing the 
second derivatives with the help of the homogeneity relation is essential to obtain a 



unique stationary solution which cannot be obtained solely from (2.30) 



We now have formally identified the quantities which determine the shape of the single- 
point velocity PDF, however, the explicit functional form of these quantities is unknown. 
Further input from the numerical or experimental side is needed to specify the PDF, 
which will be presented below, after discussing some more theoretical points. Note that 
the equations and relations derived so far are exact as no approximations have been 
made. 

2.4. Constraints on the Functional Form of the Conditional Averages 

Now that the theoretical framework is set up, we are faced with the closure problem 
of turbulence in terms of the unknown conditional averages A, /i. A, 11 and $. In this 
section, we present constraints on the functional form of these quantities, which follow 
from elementary statistical relations and statistical isotropy. These constraints are useful 
to narrow down the possible functional forms of the conditional averages. 

We start out with the pressure gradient. As already mentioned, its averaged projection 
on the velocity vector has to vanish because of incompressibility and homogeneity. This 
in turn gives an integral constraint on n(i; m 







Vp) 



dv (u ■ Vp f ) f{i 



dv vT\(v)f{v). 



(2.36) 



f We omit the i-dependence to simplify the presentation. However, we do not assume sta- 
tionarity unless explicitly mentioned. 
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With an analogous argumentation one also finds constraints for A, A and /i, 

/>oo 

- (e) = (j/M • Alt) = / dvvA{v)f{v), 



(Tr(D)> 



dv [X{v) + 2^liv)]f{v), 



and in the stationary case also for the forcing 



(e) = (u-F) = I dvv<^{v)f{v). 
la 



(2.37) 
(2.38) 

(2.39) 



A second class of constraints can be deduced from statistical isotropy, which states 
that the conditional averages are invariant under rotations e.g. 



{Vp\Yiv)=Yl{Vp\v), ReS0(3). 



(2.40) 



From this relation, it follows that the constant vector (Vp v = O) is invariant under 
rotations. However, this is possible only for the zero vector, thus 



{Vp\v = 0) =0 



n(0) = limt)- {-Vp\v) =0. 



(2.41) 



Note that the limit limt,_j.o ^ is not defined; however, the boundedness of v is enough for 
the above result. In summary, we find 



{Vp\v = O) =0 
{vl\u\v = 0) =0 

{F\v^Q) =0 



n(o) = 

A(0) = 
$(0) = 0. 



(2.42) 
(2.43) 
(2.44) 



This type of argument can also be applied to tensor-valued isotropic functions. In three 
dimensions the only constant isotropic tensor of rank two is 5ij , which leads to constraints 
on the eigenvalues 

D,j{v^0)^5,, ^ A(0) = Ai(0). (2.45) 

We can go even further and use this type of argument to obtain constraints on the 
derivatives at v = 0. For this, first note that the isotropic forms of vector- and tensor- 



valued functions (2.141 and (2.151 are also invariant under reflections. Therefore, all the 



isotropic functions considered here are actually invariant under transformations in 0(3). 
This means that in our case, the assumption of invariance under SO (3) yields invariance 
under 0(3). Now, it can be easily checked that if a tensor-valued function is invariant 
under some transformations, then so are also its derivatives. Therefore, the n-th derivative 
of D{v) at D = is a constant (n-|-2)-th-order tensor, invariant under 0(3). For n being 
odd, the derivative has to vanish, because the only constant tensor of odd order which 
is invariant under reflections is the zero tensor and it follows that 



d 

dvi. 



d 



-A 



= 



d"A 

d-u" 



d> 
dw" 



= 0, 



(2.46) 



where n is odd. With the very same argumentation applied to vectors (i.e. tensors of 
rank one) one also finds 



d"n 



dv" 



v=0 



d"A 
du" 



d"$ 



•u=0 



dw*^ 



0, 



(2.47) 



v=0 



where n is even. 

This shows that in isotropic turbulence, the series expansion of A and /i contains only 
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even powers, whereas that of H, A and <& contains only odd. These constraints have also 
to be kept in mind, when approximating the above functions with polynomials. 

2.5. A Simple Analytical Closure 
We now try to find the most simple functional forms of the conditional averages which 
still fulfil the constraints described in the previous section. By 'most simple' we mean the 



lowest-order polynomial. For A, a constant is not allowed because of (2.37) and (2.43). 



We choose a linear function and determine the pre- factor through (2.37). We obtain 



K{v) 



3ct2 



(2.48) 



where a = •y/(M^)/3 is the standard deviation of the velocity. In the stationary case, we 
further find 4>(u) — {e)/{3a'^)v. For A and /i, the most simple choice which is consistent 



with (2.38) and (2.45) is a constant 



X{v) = ^J.{v) 



Inserting this ansatz for A, fi, A into the homogeneous solution (2.33) yields 



f{v;t) 



TT a{ty 



exp 



1 



2a{ty 



(2.49) 



(2.50) 



which is the angle-integrated Gaussian distribution, i.e. the velocity vector is distributed 
according to a Gaussian 

„2 



fiv;t) 



1 



(27ra(i)2)3/2 



exp 



1 v' 



2a{ty 



(2.51) 



with standard deviation cr(t). This solution is valid for decaying as well as stationary 
turbulence, where a{t) is a monotonously decaying function in the first and a constant in 
the latter case. The above solution rests upon the homogeneity relation (2.31), which in 



the stationary case is equivalent to the evolution equation (2.32). However, for decaying 



turbulence, this evolution equation gives another possibility to obtain the PDF which 
we demonstrate in the following. For this, we choose the most simple functional form for 



the pressure gradient which is consistent with (2.36) and (2.42), namely n(w) = 0. This 
ansatz can be also viewed as an approximation which neglects the statistical dependence 
of the pressure gradient and the velocity. Furthermore, A and /x are chosen as above, 
representing statistical independence of the dissipation tensor and the velocity. With 



these assumptions, the evolution equation (2.10) simplifies to 

d 



dt 



f = - 



{e){t) 



A„/. 



(2.52) 



This is exactly the equation derived in (Hosokawa 2008) by neglecting pressure contrl 



butions and approximating the coupling of the dissipative term to the two-point PDF. 
By a change of variables ( Hosokawa||2008 ) , 

Ek^n{to)- I At' {e){t') 



1 



'(t) = ^Eu^n{t) 



1 



= r^^y 



(2.53) 



the equation takes the form of the heat equation 

d 



dr 



f = A„/. 



(2.54) 
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N' 


i?A 


L 


T Urms V 


(e) hnaxr) 


512^ 


112 


1.55 


2.86 0.543 10-=* 


0.103 2.03 



Table 1. Major simulation parameters. Number of grid points TV"^, Reynolds number based on 
the Taylor micro-scale R\, integral length scale L, large-eddy turnover time T, root-mean-square 
velocity Urms, kinematic viscosity v, average rate of energy dissipation (e). Here kmaxrj charac- 
terizes the spatial resolution of the smallest scales, where i] is the Kolmogorov length scale and 
kmax = 0.8N/2 is the highest resolved wavenumber (the factor 0.8 is due to the aliasing filter). 



The change of sign indicates that the process runs backwards in time. As stated before, 
when the fluid has come to rest, we find f{v) = 5{v), which now may serve as an initial 
condition. The solution reads as 

i.e. we have a Gaussian PDF, whose time evolution is solely determined by the mean 



rate of energy dissipation. Note that this solution is consistent with (2.51) as T{t) = 
cr(t)^/2. Equation (2.53) highlights the connection between {e){t) and (j{t) for decaying 



turbulence; because the energy dissipation reduces the kinetic energy of the flow, the 
PDF becomes narrower and narrower as a function of time. 

The results presented in this section should not be understood as a justification for 
the Gaussianity of the velocity as they have been obtained by using purely mathematical 
arguments, which do not account for any physical properties of the velocity field. For 
example, the argumentation using the homogeneous solution (presented at the beginning 
of this section) can also be applied to the vorticity field, which is known to exhibit a 
highly non-Gaussian distribution. The validity of the above results depends on how well 
the simple choices for the conditional averages agree with their functional forms measured 
in turbulence. For this comparison, we refer to the next sections, where the numerical 
results will be presented. 

With the results of this section, the Gaussian distribution can be viewed as the most 
simple solution of the equations describing the velocity PDF, which is consistent with 
the constraints presented above. In a sense it is a zeroth-order approximation to the real 
velocity PDF. Deviations of the conditional averages from their most simple forms lead 
to deviations from the Gaussian distribution. In the following sections, we will see that 
this zeroth-order approximation roughly describes the velocity PDF, but completely fails 
for the vorticity PDF. 

3. DNS results 

3.1. Some comments on the numerical simulations 

Before presenting the numerical results, some comments on the numerical scheme are in 
order. The presented results are generated with a standard, dealiased Fourier-pseudo- 
spectral code ( Canuto et aL|[T987 Hou fc Li||2007 ) for the vorticity equation. The inte- 



gration domain is a triply periodic box of box length 27r. The time-stepping scheme is 
a third-order Runge-Kutta scheme (Shu &: Osher|[l988 ) . For the statistically stationary 



simulations, a large-scale forcing has to be applied to the flow. Here, care has to be taken 
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in order to fulfil the statistical symmetries we make use of in our theoretical framework. 
After numerous tests, we chose a large-scale forcing which conserves the energy of the 
flow. This forcing has been found to deliver accurate results concerning the statistical 
symmetries. However, the results presented here do not depend on the exact form of the 
forcing, which has been verified using other forcing methods (e.g. the common method 
which holds the amplitudes of the Fourier coefficients in a band constant) . 

The numerical results presented in the following are obtained from a simulation whose 
simulation parameters are listed in tablelT] Care has been taken to produce a well-resolved 
simulation with an integral length scale which is not too large compared with the box 
length of 27r. We have found that the statistical symmetries, especially isotropy, are not 
valid for a snapshot of the velocity field at a single point in time, but can be obtained 
through time- or ensemble-averaging, respectively. This is due to long-range correlation of 
the velocity and is not as pronounced for the vorticity, a short-range correlated quantity. 
To ensure the validity of the statistical symmetries, the averages for the stationary case 
have been obtained through spatial and temporal averaging over more than 150 large- 
eddy turnover times, whereas for decaying turbulence a spatial and an ensemble average 
over 12 realizations has been applied. 

It may be doubted that the results obtained from a single numerical set-up are of 
universal nature, for example there may be influences of the imposed (periodic) boundary 
conditions. This critique, of course, also applies to any experimental result, so that a 
compilation of data from many different numerical and experimental set-ups would be 
useful to discuss the issue on a more general level. However, the results presented in the 
following will make a good point that there are deviations from the Gaussian shape in 
both stationary and decaying turbulence. 

3.2. Stationary Turbulence 



The homogeneous and stationary solutions (2.331 and (2.351 of the PDF equations allow 



a detailed analysis of the connection between the functional shape of the PDF and the 
statistical correlations of the different dynamical quantities, represented by the condi- 
tional averages H, A, $, A and fi. This analysis is presented in this section for the case 
of stationary turbulence. 

To start with, we consider the diffusive term, shown in figure [ll It exhibits a negative 
correlation with the velocity, indicating that a fluid element will on average be decelerated 
due to viscous forces. The pressure term, however, displays an interesting zero-crossing. 
This means that a fluid particle will on average be decelerated for low values of velocity 
and accelerated for high values of velocity. If the pressure term is non-vanishing, it has 



to exhibit this zero-crossing, because the integral constraint (2.361 can be fulfilled only 
through a zero-crossing or vanishing of H(w). Our numerical investigations show that the 
pressure contributions to the single-point velocity PDF cannot be neglected, in contrast 



to the assumptions of some recent theories ( Tatsumi fc Yoshimura|2004 Hosokawa|2008 ). 



This is especially true for high velocities, where we see a strong ii-dependence of 11 and 
hence a clear statistical dependence of the pressure gradient and the absolute value of 
the velocity. 

The conditionally averaged forcing term in figure [l] was computed based on the con- 



ditional balance (2.34 1. This is because the forcing in our numerical scheme is implicitly 
applied to the vorticity field and is not available as an additive field that we can aver- 
age conditionally. However, some careful tests have been performed to ensure that this 
balance actually holds. An interesting feature of $ is that it seems to saturate for large 
values of v in contrast to H and A, suggesting that the forcing plays a minor role for 



the statistics of high velocities. In (Falkovich fc Lebedev||1997) it has been argued that 
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0.5 1 1.5 



Figure 1. The conditional averages Ti-{v), A{v), ${v) and the PDF f{v). The w-axis is scaled 
with the standard deviation of the velocity a = Urms = \/{ii^/3). The conditional averages 
show an explicit ii-dependence indicating statistical correlations. 




Figure 2. The conditional averages /i(i'), A(t;), (e|u), (i^uj^\v) and the PDF f{v). The 
conditional averages show an explicit u-dependence indicating statistical correlations. 



the tails of the velocity PDF are related to those of the forcing. Our numerical results, 
however, suggest that the main contribution conies from the dynamical effects of the 
pressure gradient and the velocity diffusion. Furthermore, the saturation of $ suggests a 
simple relation between 11 and A for high velocities, namely an approximately constant 
offset. This shows that for high velocities, the statistical dependences of the pressure gra- 
dient and the velocity diffusion on the velocity arc strongly related. A natural question 
concerning the conditionally averaged forcing is whether its form depends on the type 
of forcing. We have tested two other common forcing methods (these methods manip- 
ulate the vorticity in a Fourier band by holding either the Fourier coefhcients or their 
amplitudes constant), finding that the general shape of $ does not depend on the forcing 
method. 

Now let us proceed to investigate the conditional energy dissipation tensor, its eigen- 
values are shown in figure I2I together with the PDF f{v). Both eigenvalues have a similar 
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Figure 3. Comparison of the PDF f{v) with the evaluation of (2.331. The reconstructed PDF 
collapses with the directly estimated PDF; therefore, the latter is marked by black dots to 
indicate its presence. The excellent agreement demonstrates the consistency of the theoretical 
results. A comparison with an angle-integrated Gaussian shows that, although being close to 
Gaussian, significant deviations exist. 



dependence on v, they are almost constant for small and medium velocities, but display a 
strong ti-dependence for high velocities. This shows that the occurrence of high velocities 
is statistically correlated with strong dissipation and enstrophy. 

The conditional averages shown in figures [l] and [2] have been tested to obey the con- 
straints presented in § 2.4 The validity of the constraints &l v — Q can be seen in the 



figures, whereas the integral constraints have been verified numerically. As a test for 



consistency, the conditional averages are inserted into (2.331 to calculate the homoge- 



neous PDF, and the result is presented in figure[3j An excellent agreement with the PDF 
directly estimated from the numerical data is found. The same results can be obtained 



using (2.35) as it is equivalent to (2.33) in the stationary case. A comparison with an 



angle-integrated Gaussian (2.50) shows that, although the PDF is similar to Gaussian, 



significant deviations occur. Especially notable is the sub-Gaussian tail, which has also 
been found in experiments and other numerical simulations ( Vincent fc Meneguzzi|199r 
Noullez et a/.|[l997i[Gotoh et al.pm^ 



It is well known for a long time that the velocity PDF is similar to a Gaussian distribu- 
tion. However, this fact (coming from experiments and simulations) is a priori not clear 
at all. One can certainly arrive at a Gaussian distribution by neglecting certain terms or 
assuming the most simple forms for the conditional averages (see § 2.5). However, as can 



be clearly seen in figures [T] and [21 those assumptions are violated. Furthermore, some of 
the arguments can also be applied to the vorticity which is known to exhibit a highly 
non-Gaussian PDF. This shows that one has to be very careful with the utilization of 
such simple arguments. Therefore, here we rather use the measured conditional aver- 
ages to understand the similarity of the velocity PDF to the Gaussian distribution. The 
question sought after is, how it is possible to find strong statistical correlations for the 
various conditional averages and the velocity, but still only moderate deviations from 
Gaussianity for the PDF. 

As already stated, there are clear deviations of the conditional averages from the simple 
forms suggested in § |2.5| We see that for small and medium velocities, the assumption 
of linear A (or equivalently of linear $ and vanishing H) and constant A, // can be seen 
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as a rough approximation. Because of this, the PDF stays close to a Gaussian in this 
range of velocities. However, for high velocities we see a strong w-dependence of A, fi and 
a functional form of A which is clearly not linear. Still, only moderate deviations from a 
Gaussian tail occur. To understand this, let us examine again the relation between the 
PDF and the conditional averages: 

j\f r ,Hv') + ^ti'W) 

We see that only the quotients A/A and /^/A enter the exponential function. As —A and 
A both are bended upwards for high velocities, their quotient might still be linear, i.e. 
the deviations from the simple forms of § |2.5| might compensate each other. Indeed, this 
is roughly true as can be seen in figure [TT] which also shows that A and ^ are roughly 
equal. Furthermore, a comparison with the case of vorticity is shown, which will be 
discussed later. Assuming for the moment a linear A/A and jjl/X = 1 yields a Gaussian 



shape for the exponential factor in (3.1 1. In this case, the second factor, 1/A, is then 
responsible for (slight) modifications of the Gaussian shape. Therefore, it appears that 
the compensation of A and A to yield an approximately linear quotient A/A is crucial 
for the similarity of the velocity PDF to the Gaussian distribution. For example, if one 



artificially sets X — ^ = (e)/3 as in § 2.5 but still uses the measured and highly nonlinear 
A for the reconstruction, one obtains a PDF which decays much faster than a Gaussian. 
To discuss the issue of the sub-Gaussian tails, we need to take into account the de- 
viations of the quotients from a linear and constant function. A detailed analysis shows 
that the deviations from Gaussianity cannot be solely attributed to the pre-factor 1/A, 
but are also due to the aforementioned deviations of the quotients. What makes the issue 
even more complicated is that the relative amount of contribution to the deviations by 



the two factors in (3.1) seems to vary with the Reynolds number, as has been checked 
with a number of simulations not presented here. This shows that the detailed properties 
of the PDF of the velocity are due to a subtle interplay of the statistical dependences 
of dynamical quantities such as the dissipation, the velocity diffusion or the pressure 
gradient. 

3.2.1. Approximating the Dissipation Tensor 

In view of the complicated statistical dependences, it is natural to ask for approxi- 
mations which simplify the matter but still describe the velocity PDF reasonably well 
e.g. maintain the sub-Gaussianity of the velocity PDF. It turns out that the conditional 
dissipation tensor allows simplifications which contain some interesting physical insights. 
It appears from figure [2] that the eigenvalues of D have a similar functional form, sug- 
gesting the approximation X{v) « fJ-iv). Thus, all eigenvalues of D are approximately 
equal, which basically states that the conditional dissipation tensor does not contain any 
directional information and is proportional to the identity matrix 

jTr(D)..,4. 

In this approximation, D depends only on the magnitude of the velocity, i.e. the dissi- 
pation tensor i^AA and the direction of the velocity vector are uncorrelated. Note that 
this property is not a consequence of the statistical isotropy, but rather a special prop- 
erty of the velocity in turbulent flows (this will become more clear when we examine the 
dissipation tensor of the vorticity). 

The equality of the eigenvalues may be motivated with a simple scale-separation argu- 
ment. The velocity itself is known to be a rather smooth vector field varying quite slowly 



D,,{v) « -Tr(D) S,j = -{e + iyuj^\v) 6,j. (3.2) 
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in space. Compared to that, the gradients vary on a much shorter scale (see figure [9]). 
The numerical results suggest that quantities like, for example, the dissipation depend 
on the absolute value v. If we now assume that the dependence on v is crucial, but due 
to scale separation no correlation between the direction of u and the velocity gradient 
tensor exists, we may consider a projection on m as a projection on a random direction. 
This immediately gives 

{v{k^uf\v) = {vukk^u\v) = l{e + vuj''\v), (3.3) 

6 

i.e. averaging the projection of a tensor on a random direction yields one-third of the 
conditionally averaged trace of the tensor. This in turn yields \{v) = m(^)j which is 
suggested by the numerical results as an approximation. 

As the trace is the central quantity, one may wonder if further simplifications are pos- 
sible. To this end, we consider the average proportionality of enstrophy and dissipation in 
homogeneous turbulence, (e) = {vuP'), which comes due to (Ap) = 0. It is now straight- 
forward to ask whether this balance also holds in a stronger sense, when the conditional 
average is considered. The numerical evaluation of these terms is also presented in figure 
|2J supporting that this stronger equality actually holds. As a consequence, the approxi- 
mation \{£ + viJ^\v) ss {e\v) seems reasonable. Taking this all together, the conditional 
dissipation tensor can be approximated as 



D.,,{v)^-{e\v)5,,, (3.4) 



1 

3^ 
depending only on the conditionally averaged rate of kinetic energy dissipation. 

It may be hypothesized that the approximate relations so far only hold asymptotically 
with e.g. increasing Reynolds number. By comparing simulations with different Reynolds 
numbers (up to Re\ « 200), we were not able to clearly verify (or falsify) this statement; 
hence, simulations with significantly higher Reynolds numbers are needed to answer this 
question. 

In the light of these simplifications, it is interesting to reconsider the functional form 



of the velocity PDF. Equations (2.161 and (3.1) together with X{v) « ijl{v) ~ (e|i;)/3 
instantly yield 

M r , ,„ (iyu-Au\v') 



f{v) ^ i-TT exp / dv' 3 



[e v 



N [\ ,^ {u-Vp~u-F\v') 

■exp / dv 3 , . , , (>J-5) 



{e\v) "/„„ {eW) 

indicating that the shape of the velocity PDF may be solely determined by the velocity 
diffusion and the rate of energy dissipation. As the diffusive term balances the sum of 
the pressure and the forcing term, the PDF is equivalently determined by how energy is 
injected into the system, redistributed by non-local pressure effects and finally dissipated. 
Figure |4] compares this approximation with the measured PDF and a Gaussian. We see 
that the above formula for the velocity PDF holds quite well for the core of the PDF, but 
leads to deviations in the tail. This (in addition to figure [2]) shows that the assumption 
of absent correlation between the dissipation tensor and the direction of the velocity 
becomes less valid for high velocities. Nevertheless, this approximation describes the 
velocity significantly better than the Gaussian distribution. 

We note here that the closeness of the velocity PDF to a Gaussian can be attributed 



On the velocity distribution in homogeneous isotropic turbulence 



17 




Figure 4. Comparison of the PDF f{v) wit h th e angle-integrated PDF of the approximation 

(3.51. 




Figure 5. Velocity PDF f{v;t) for decaying turbulence at time instants separated by 0.35T, T 
being the large-eddy turnover time of the initial condition. The insets show the PDFs normalized 
to unit variance, indicating the existence of a self-similar regime. 



to the approximate validity of 



{lyii ■ Au\v) 
{s\v) 



V 



(3.6) 



It describes the relation of the dissipation and the velocity diffusion, which is actually the 
dynamical effect responsible for the dissipation. It is the theoretical framework presented 



in the previous sections and in particular the formula (3.1) which allows to relate the 



observation of closely Gaussian velocity PDF to the non-obvious relation above. 

3.3. Decaying Turbulence 

We now study the shape and evolution of the velocity PDF in the case of decaying 
turbulence. The motivation for this is two fold. First, it has recently been hypothe- 



sized by Hosokawa (2008) that decaying turbulence might be profoundly different from 
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10" 




Figure 6. Comparison of f{v) (directly estimated and reconstructed from ( 2.33 1) together with 
an angle- integrated Gaussian. A deviation from the Gaussian shape is also found in the decaying 
case. The structure of the conditional averages is similar to the case of stationary turbulence. 
The averages and the PDFs are computed from an ensemble for a decay time of lAT. 




Figure 7. Sum of the p ressure and diffusive terms for decaying turbulence together with the 
analytical relation ( |3.13[ ) in straight lines. Slight deviations are visible after O.TT, which are 
negligible later on. 



forced, stationary turbulence. As mentioned in § 2.5 Hosokawa found Gaussian solutions 
for decaying turbulence by approximate closure assumptions. Second, in the work by 
Falkovich & Lebedev (1997), the argument for sub-Gaussian tails of the velocity PDF 



relies strongly on the statistics of the external forcing. This effect, of course, is absent 
for decaying turbulence. 

Figure p^ shows the velocity PDF for decaying turbulence f{v;t) for different points 
in time. The results were obtained by averaging over an ensemble of 12 independent 
simulations to increase the statistical quality. The variance of the PDFs decreases as a 
function of time as the kinetic energy of the system is dissipated. However, when rescaled 
to unit variance, the PDFs collapse, indicating a self-similar regime. Hence, once the PDF 
is specified at the beginning of this regime, it is determined subsequently according to 



a{t) \ cr{t) 



(3.7) 



where cr(i) may be related to the kinetic energy or the rate of energy dissipation according 
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Figure 8. Kinetic energy for decaying turbulence as a function of time. After a transient the 
energy decays algebraically with an exponent close to 1.4. The inset shows the rate of energy 
dissipation, which also exhibits an approximately algebraic decay with an exponent close to 2.4. 



to (2.531. Turning to a closer investigation of the functional form of the PDF in this 



regime still reveals deviations from Gaussianity, as can be seen in figure [6j Here, also the 
conditional averages of the pressure term, the diffusive term and the terms related to the 
conditional dissipation tensor are shown. They all show a functional form coinciding with 
the observations in stationary turbulence, indicating that the statistical correlations and 
the shape of the PDF do not fundamentally differ for decaying turbulence. As a check for 



consistency, the PDF is computed from the conditional averages by (2.331, which again 
performs very well. On the basis of these results, we can conclude that sub-Gaussian 
velocity PDFs are also found for decaying turbulence, which can be tracked down to the 
interplay of statistical correlations already discussed in the stationary case. 

It is interesting to study this self-similar decay more deeply in the framework of the 
PDF equation (2.30). For decaying turbulence, it is particularly useful to employ the 
method of characteristics. Let V{t, vq) denote a characteristic curve which starts from 
1^0, i.e. V(tQ, Vq) — Vq. The method of characteristics then yields the following equations: 



di 



Vit,vo) 



piV(t,VQ) 



t) 



mv,t)+Aiv,t)i^y^,^^^^, 

d_ 
dv 



(n(«,i) + A(t-,t)) 



I{V{t,vo)-t), 



(3.8) 
(3.9) 



v=V(t,va) 



of which the latter is easily integrated. We obtain the evolution of / along the charac- 
teristic curves: 



f{V{t,VQ)]t) = f{vo]to) exp 



- / di' 

to 



d_ 

dv 



'n{v,t') + K{v,t')) 



v=V{t',vo) 



(3.10) 



Of course, we are more interested in the temporal evolution of /(w; t) instead of /(V^(t, VQ)\t). 
This mapping can be achieved with the inverse function of y(i,wo), which is defined by 
V^^{t,V{t,VQ)) = Vq. Thus, the temporal evolution of the PDF of the magnitude of 
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velocity is given by f{v;t) — f {V {t,V~^ {t, v)); t) , such that we obtain 



/(z;;i) = /(y-i(t,v);io)exp 



dt' 



to 



_d_ 



-(nK,i') + AK,i')) 

Jt,'=y(t'.\/-i(t.D)) 

(3.11) 
The interesting fact about this result now is that we obtain self-similar solutions of the 



form (3.7) if and only if the characteristic curves are of the form 



V{t,vo) 



(Jito) 



Vo, 



which due to (3.8) gives rise to the relation 



U(v,t)+A(v,t) =u — In 
dt 






1 {e){t) 

2 Ekin{t) 



V. 



(3.12) 



(3.13) 



This means that we only obtain self-similar solutions if the sum of the conditional averages 
related to the pressure gradient and the Laplacian of the velocity is a linear function 
of V, with a negative slope proportional to the ratio of the rate of energy dissipation 
and the kinetic energy. Note that this non-trivial relation is a direct consequence of the 
observation of a self-similar regime of the velocity PDF, no matter whether it is Gaussian 
or not. The relation is also checked with our numerical results with some examples shown 
in figure ItI Apart from the beginning of the decay phase, the analytical relation (3.13) 
is in very good agreement with measured conditional averages. 

To characterize the statistics of decaying turbulence further, we study the temporal 
evolution of the kinetic energy and the rate of energy dissipation. The decay rate of 
energy is a central question in turbulence research, which still remains a point of debate. 



We refer the reader to Karman & Howarth (1938); Kolmogorov (1941); Batchelor & 



Townsend (1948); Saffman (1967); George (1992) for a detailed account of the theoretical 



aspects of this question as well as some classical experimental results. While a power-law 
decay of the kinetic energy is widely accepted, the prediction for the precise numeric 
value of the decay exponent varies. Detailed numerical investigations on this have been 
presented recently in (Ishida et al. 2006 Perot 2011). Apart from the dependence of 



the decay exponent on some fundamental theoretical issues concerning the Loitsyansky 
integral, it has to be assumed that it depends on properties of the initial conditions 
( George 1992 Ishida et al. 2006 ) . As we have not investigated a larger class of initial 



conditions, the following results shall characterize the simulations rather than making a 
point for or against some of the cited theories. 

The kinetic energy for a single run is shown in figure [8J After a short transient period 
of about IT (T is the large-eddy turnover time of the initial condition) , the kinetic energy 
decays algebraically with an exponent close to 1.4, i.e. Ekin{t) ^ t~-^-'^, which is in good 
agreement with the results presented in ( Ishida et a^.||2006 ) and close to the Kolmogorov 
prediction of Ekinit) ~ t~^'^/^. Accordingly, the kinetic energy dissipation decays with 
an exponent close to 2.4 in this regime, as shown in the inset of figure [8j The self-similar 
decay together with an algebraic evolution of the variance indicates a particularly simple 
evolution in time. It is also interesting to note that the slight deviations from linearity of 
the conditional acceleration observed in figure [7] lie within the short initial transient. It 
seems as if the turbulence needs a short relaxation time to switch from the forced regime 
to the decaying regime. 
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Figure 9. Volume rendering of the magnitude of the velocity field (left top), the trace of the 
dissipation tensor j^Tr(AA^) = ^{e + vlj^) (also known as pseudo-dissipation, left bottom), the 
magnitude of the vorticity field (right top) and the trace of the enstrophy dissipation tensor 
!/Tr [Vaj(Vu))^] (right bottom). The colours blue, green, yellow, red in that order indicate 
increasing amplitude. The vorticity and both dissipation tensors tend to organize into entan- 
gled small-scale structures forming a complicated global structure. The velocity appears more 
unstructured but displays long-range correlations. The visualizations have been produced with 
VAPOR (www.vapor.ucar.edu). 



3.4. Comparison with the Vorticity 

To demonstrate that some of the results are unique features of the velocity statistics, we 
draw a comparison with the statistics of the turbulent vorticity, which has been studied 
more extensively in fWilcz ek &: Friedrich||200 9) . A first glance at snapshots of the fields 
reveals a fundamental difference between velocity fields and vorticity fields as shown in 
figure [9J While the vorticity tends to be organized in small filamentary structures, the 
velocity seems to be more structureless and smeared out. Recall that the velocity may 
be calculated from the vorticity with the help of the Biot-Savart law, which may be 
interpreted as some kind of a smoothing filter. On the level of the single-point statistics. 
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the most striking difference is that the vorticity displays a highly non-Gaussian PDF 
with pronounced tails, which is closely related to the coherent structures present in the 
flow. Moreover, as for example reasoned in ( Tennekes fc Lumley|1983 ), in the case of the 
turbulent vorticity the average enstrophy production and enstrophy dissipation tend to 
cancel, which means that the external forcing has a negligible effect on the statistics of 
the vorticity at high Reynolds numbers. 

The evolution equation for the vorticity PDF /o(ri;a;,i) can be derived analogously 
to the PDF equation for the velocity. For homogeneous turbulence, it takes the form 

^/n - -Vo • [(Su; + uAL^ + VxF\VL)fn]. (3.14) 

Exploiting homogeneity, the diffusive term may also be rearranged, such that the PDF 
equation takes the form 



d_ 
di 



In = -Vs7 • [{SoJ + V X F |f2) /o] - 
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(3.15) 



d^idflj \ dxk dxk 

Hereby the conditional enstrophy dissipation tensor enters. Because of statistical isotropy 
the quantities in these equations take the form 



(VxF|rj) = $n(n)fi, 



Dn■.^J{n) 



2 

duji dujj 
dxk dxk 



n 






^^m = {u ■ {V X F)\n) , 



n^ 



(3.16) 

(3.17) 

(3.18) 
(3.19) 



where Jl and a> are the unit vectors belonging to $7 and lj. The rate-of-strain tensor 
may, due to incompressibility, be characterized by a single scalar function I](f2) related 
to the enstrophy production. The conditionally averaged vortex stretching term can then 
be expressed as 

(Sa;|n) = (S|0)0 = S(^)r20. (3.20) 

The eigenvalues of D^ can be computed using relations analogous to (2.261 and (2.27). 
The isotropic forms of the evolution equations and the homogeneity relation (analogous 
to ([230|-([232l)) are 



d ~ 


d 


0== 


d 


d ~ 


d 
90 



(EO 
An 






S0 + $ 



o 



$0) k 



ao2 



fa 



Ao/n 



902 



An/n- 



(3.21) 
(3.22) 

(3.23) 

A„(n',*) ■ ''■^'" 

The reconstruction of the PDF using this formula again works perfectly well, as shown in 
figure 10 For stationary turbulence, the probability current in (3.21) has to vanish and 
therefore also the conditional average in (3.14). However, it turns out that the effect of 



Solving (3.22) yields a relation between the PDF and the conditional averages 
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Figure 10. The PDF of the vorticity and the conditional averages which determine its shape. 
Note that the forcing term vanishes and therefore An and E Q cancel. The abscissa is normalized 
by the standard deviation of the vorticity, whereas the ordinate is normalized using the average 

of en ~ i'/2Tr ( [Voj + (Vc<j)^l ) and the standard deviation of the vorticity an. Note that not 



An but An/4 is shown in the right figure. 



the external forcing may be neglected here, as can be seen from figure \T0\ Therefore, the 
conditionahy averaged vortex stretching term and the diffusive term tend to cancel almost 
identically. This means that the balance between enstrophy production and dissipation 
holds even under conditional averaging, which has already been shown in ( Novikov|1993 
Novikov & Dommermuth 1994 Wilczek & Friedrich 2009). As a consequence, we may 



regard the enstrophy production as an internal process, which is decoupled from the 
external forcing. This interpretation is consistent with the common assumption that the 
small-scale properties of turbulence become independent of the forcing mechanism given 
a sufficiently high Reynolds number. As demonstrated in the preceding paragraphs, the 
forcing term may not be neglected when determining the stationary PDF of the velocity. 
The conditional averages are shown in figure [lO] together with the vorticity PDF. 
They display a strong dependence on 51, especially the eigenvalues which do not show 
a plateau for small values of fl like in the case of velocity. An important observation 
here is that, in contrast to the velocity statistics Ao 56 /ijj, but rather An sa 4/isi- This 
shows that the tensor Do depends on the direction of Jl, i.e. the enstrophy dissipation 
tensor v^u (Va;) and the direction of the vorticity Jl are obviously correlated, in 
contrast to an almost vanishing correlation in the case of velocity. One may hypothesize 
that the simple scale-separation argument, which works in the case of the conditional 
dissipation tensor, breaks down for the conditional enstrophy dissipation tensor, because 
there is no clear scale separation between the vorticity and its gradients. This can be 
observed in figure |9] where a snapshot of the velocity, vorticity and the traces of the 
corresponding dissipation tensors are shown. It appears that the fields of the vorticity 
and both dissipation tensors consist of structures of similar size, whereas the velocity 
field consists of fairly large patches. Furthermore, the structureflf] of vorticity and the 
dissipation tensors are oriented in almost the same direction, but no obvious correlation 
of the orientation of the structures of the velocity and the other fields is visible. However, 
the amplitudes of the fields are correlated. Therefore it is understandable that D and Dq 
both depend on the magnitude of the corresponding variable, but the dependence on the 

f The vorticity field consists of tube-like structures, whereas the pseudo-dissipation field 
!/Tr(AA''") — |(e + ulu'^) consists of structures, which are a superposition of sheet-like and 
tube-like structures. 
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Figure 11. The quotients A/ A and ^/A appearing in (2.331 and (3.241 for the case of velocity 

(left) and vorticity (right). 



direction of the variable {v and 17 respectively) is weak in the case of velocity and strong 
in the case of vorticity. 

We have already noted that the arguments in 



2.5 which lead to a Gaussian PDF, 



could, in principle, also be applied to the vorticity. Figure 10 shows a comparison of fo, 
with an angle-integrated Gaussian. This exemplifies that such types of arguments can 
fail completely and therefore cannot be used as proofs of Gaussianity. To understand 
the reasons for a strongly non-Gaussian vorticity PDF we again examine the quotients 
Ao/Ao and ^o/An, depicted in figure 11 We see that Ao/Ao is clearly not hnear and has 



a complicated dependence on fi, showing that there is a complicated relation between 
vorticity diffusion and the enstrophy dissipation tensor. The quotient of the eigenvalues 
also shows a dependence on il, corresponding to the fact that jii^ o^ \q_. K detailed 
analysis shows that the complicated form of the vorticity PDF depends crucially on the 
non-trivial form of both quotients, i.e. it cannot be attributed to one of them. Comparing 
the quotients for the case of velocity and vorticity, we see that there is much more non- 
trivial behaviour in the vorticity case. This leads to a highly non-Gaussian vorticity PDF, 
whereas the velocity PDF displays only moderate, but still significant, deviations from a 
Gaussian distribution. 

In this section, we have examined the difference between the velocity and vorticity 
PDFs, using the conditional averages which determine these PDFs. We have seen that 
the approximate relations A ~ /i and A/A ^ v, which have been found for the velocity, 
are clearly violated in the case of vorticity. Overall, we can say that the conditional 
averages or the statistical dependences which determine the shape of the vorticity PDF 
are significantly more complicated for the vorticity and hence yield a highly non-Gaussian 
PDF. 



4. Summary and Conclusions 

To sum up, we studied the statistics of the single-point velocity PDF within the frame- 
work of the LMN hierarchy. Combined with conditional averaging and the use of sta- 
tistical symmetries, this framework provides a clear-cut identification of the quantities 
which determine the details of the velocity PDF. The theory identifies the conditional 
diffusion of velocity and the conditional dissipation tensor in the homogeneous case as 
the central quantities of interest. In the case of stationary turbulence, the diffusive term 
is equivalently represented by the conditional pressure gradient and the external forcing. 
Exact expressions for the velocity PDF are presented in terms of these functions, demon- 
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strating that the closure problem of turbulence may be expressed in terms of a priori 
unknown statistical correlations of different dynamical contributions directly related to 
the Navier-Stokes equation. 

In an analytical treatment, we suggest simple closure approximations that comply 
with the functional constraints on the conditional averages leading to Gaussian PDFs 
for both stationary and decaying turbulence. It turns out that the presented closure 



approximations contain the results recently derived in ( Hosokawa 2008 ) , explaining which 
correlations have to be neglected to yield Gaussian statistics. 

DNS simulations of stationary and decaying turbulence are then used to study devia- 
tions from Gaussianity. It turns out that the velocity PDF displays slight deviations from 
a Gaussian distribution with sub-Gaussian tails. A detailed investigation of the terms 
arising in the theoretical framework allows to investigate why only moderate deviations 
from Gaussianity occur in spite of pronounced correlations. Unlike suggested in recent 
theories, we show that the pressure contributions may not be neglected. Furthermore, the 
correlations of the velocity with the external forcing do not seem to be the main contrib- 
utor to the deviations from Gaussianity. The conditional energy dissipation tensor turns 
out to display strong correlations for high values of the velocity, which contribute to the 
deviations. The numerical results additionally suggest that this tensor is approximately 
diagonal with identical eigenvalues related to the conditional rate of energy dissipation. 
The investigation of decaying turbulence gives similar results, indicating that forced and 
decaying turbulence do not differ in a fundamental way. Interestingly, a self-similar range 
with an algebraic decay of the kinetic energy has been observed, which considerably sim- 
plifies the description of the evolution of the PDF during the decay phase and implies 
a simple relation for the conditionally averaged right-hand side of the Navier-Stokes 
equation. 

To highlight the genuine features of the velocity statistics, some comparisons to the 
vorticity PDF are drawn. The general outcome here is that the vorticity is stronger corre- 
lated with the terms arising in the vorticity equation that determine the local dynamics 
such as enstrophy production and dissipation. The external forcing has a negligible effect 
as has been found in various other studies. The conditional enstrophy dissipation ten- 
sor cannot be assumed to be diagonal, in contrast to the conditional energy dissipation 
tensor, indicating the presence of non-negligible directional correlations. This fact may 
be related to the presence of coherent structures in the case of the vorticity and their 
absence in the case of the velocity; however, a concise description of this fact is still 
lacking. 

In conclusion, the approach presented above gives a comprehensive characterization 
of the shape and evolution of the single-point velocity PDF. Although the introduction 
of conditional averages does not solve, but rather reformulates the closure problem of 
turbulence, valuable insights have been obtained highlighting the statistical correlations 
of the different dynamical influences and the velocity. Regarding a statistical theory 
of turbulence, it is interesting to see how PDFs may be expressed in terms of these 
influences. A theoretical derivation of the functional shape of these quantities remains a 
challenging task for the future. 
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